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ABSTRACT 


Through  the  introduction  of  a  new  parameter  the  radiation  pattern 
of  a  phased  array  using  digital  phase  shifters  for  beam  steering  can  be 
made  periodic  and  pattern  considerations  can  be  confined  to  a  very  narrow 
range  of  the  scan  angle.  Arrays  designed  on  the  basis  of  a  least  mean- 
squared  pattern  error  are  shown  to  require  the  smallest  absolute  values 
of  a  phase-index  function.  A  systematic  and  easy-to-apply  method  is 
developed  for  reducing  the  peak  sidelobes  of  nonreciprocal  phased  arrays 
by  one-step  phase  adjustments  in  certain  elements.  The  amount  of  possi¬ 
ble  sidelobe  reduction  depends  on  the  scan  angle  and  the  size  of  phase- 
quantization  steps.  Typical  examples  for  4-bit  phase  shifters  show  that 
reductions  in  excess  of  9  dB  are  possible  for  certain  main-beam  directions. 
By  the  appropriate  choice  of  a  weighting  function,  the  same  technique  can 
be  readily  used  to  suppress  noise  or  interference  coming  from  any  given 
direction. 
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I.  introduction 


Increased  attention  has  recently  been  directed  to  the  development 
of  phased  arrays  which  utilize  digital  phase  shifters  for  beam  steering. 
It  Is  well  known  that  a  digitally  controlled  phased  array  designed  on 
the  premise  of  continuous  phase  shifts  produces  undesirable,  high  side- 
lobes  because  of  the  existence  of  phase-quantization  errors.  [1,2]  These 
sidelobes  diminish  when  a  large  number  of  bits  are  used  In  the  digital 
phase  shifters.  However,  constraints  on  cost  and  system  complexity  can 
not  allow  arbitrarily  small  steps  In  phase  changes.  It  Is  therefore  Im¬ 
portant  to  study  the  dependence  of  an  array  pattern  on  the  quantized 
phase  errors  and  to  search  for  effective  ways  for  reducing  peak  sidelobes. 

The  purpose  of  this  report  Is  threefold.  First,  It  will  be  shown 
that  by  Introducing  a  new  parameter  In  place  of  the  scan  angle  an  array 
pattern  can  be  made  periodic  with  respect  to  the  new  parameter  Inde¬ 
pendently  of  the  element  spacing  or  the  quantized  phase  step.  This  sim¬ 
plifies  the  study  considerably  since  It  will  only  be  necessary  to  consider 
the  array  pattern  within  a  very  narrow  range  of  the  scan  angle.  Second,  1, 
will  be  demonstrated  that  the  simple  expediency  of  minimizing  the  mean- 
square  error  of  the  array-pattern  function  does  not  yield  uniformly  low 
sidelobes  throughout  the  visible  region;  some  sidelobes  may  become  In¬ 
tolerably  high.  Third,  a  new  technique  will  be  presented  for  reducing 
the  peak  sidelobes  of  digitally  controlled  radar  phased  arrays.  This 
technique  takes  advantage  of  the  nonreciprocal  nature  of  latching  ferrite 
phase  shifters.  It  Is  systematic  and  does  not  require  random  trials. 
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Examples  will  show  that  the  peak-sidelobe  level  obtained  on  the  basis  of 
a  minimum  squared  phase  error  can  be  further  reduced  by  9  or  more  dB  in 
certain  cases.  By  suitably  controlling  a  weighting  function,  this  tech¬ 
nique  can  be  readily  used  to  suppress  noise  or  interference  coming  from 
any  given  direction. 
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•ii-i — PE.RIODICITY  OF  DIGITALLY  CONTROLLED  ARRAY  PATTERNS 


Consider  a  linear  array  of  2K+1  elements  with  uniform  spacing  d, 

as  shown  in  Fig.  1.  Let  Ik  and  -27rn(k)/M  denote,  respectively,  the 

excitation  amplitude  and  the  phase  shift  in  the  kth  element,  where  n(k) 

is  an  integer  and  M*2i  for  an  i-bit  phase  shifter.  The  far  field  of  the 
array  is 


K 

E(6)  =  l  I,  exp{ j2w[--  sin  6  -  S^.]}  , 
k=-K  A  M 

(1) 

where 

of  the 

6  is  measured  from  the  broadside  direction.  If  0  is 

o 

main  beam,  (1)  can  be  rewritten  as 

the  direction 

K 

S(u)  -  I  I.  exp(jku)  exp[j  ~f(a,k)]  , 

k«-K  n 

(2) 

where 

2tt 

U  -  —  (sin  6  -  sin  0  ) 

A  O 

(3) 

Md  .  A 
a  •  sin  0 

A  o 

(4) 

f(a,k)  =  ak  -  n(k). 

(5) 

From  (2),  it  is  clear  that  the  condition  f(a,k)  =  0  yields  the  radiation 
pattern  obtained  with  continuous  phase  shifts  and  that  the  deterioration 
of  the  pattern  due  to  quantized  phase  shifts  is  uniquely  determined  by  the 
quantity  27rf(a,k)/M.  Moreover,  since  0q  is  a  single-valued  function  of  a 
in  the  region  -tt/2  <  6o  <  */2,  a,  instead  of  0q,  can  be  used  to  denote 
the  main-beam  direction.  For  convenience,  we  shall  call  f(a,k)  a 
phase-index  function. 
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Let  -2nno(a ,k)/M  be  the  phase  shift  in  the  kth  element  which  pro¬ 
duces  a  pattern  pointing  in  the  a-direction  with  a  prescribed  sidelobe 
level.  Then,  if  we  define  n  (o+m,k)  for  the  new  main-beam  direction 

i 

a+m  as 

I 

no(a-hn,k)  =  no(a,k)  +  mk  ,  (6) 

we  hrve  from  (5)  and  (6), 

fo(o+m,k)  =  fQ(a,k)  .  (7) 


Equation  (7)  implies  that  the  pattern  with  beam  direction  a+m  has  the 
same  structure  and  the  same  prescribed  sidelobe  level  as  that  with  beam 
direction  a.  In  case  a  lies  in  the  region  0  <_  a  1/2,  we  define 


no(1"a»k)  “  k  "  n0(“»k)  •  (8) 

If  the  array  has  a  symmetrical  amplitude  excitation  and  an  antisym- 
metrical  phase  distribution;  that  is,  if  Ifc  «  I_k  and  no(a,k)  -  -  no(a,-k), 
we  obtain,  on  substituting  (8)  In  (5) 


f0U-«,W  »  -  fQ(a,k) 


-  f 

o 


(a*-k). 


(9) 


From  (2.)  and  (9),  it  is  seen  that  the  pattern  with  beam  direction  1-a  is 
symmetrical  to  the  pattern  with  beam  direction  a  with  respect  to  u=Q.  As 
a  result,  the  consideration  of  the  array  radiation  pattern  can  be  restricted 
to  the  range  0  a  1/2,  from  which  the  radiation  pattern  for  any  arbitrary 
value  of  a  can  be  derived.  The  corresponding  range  in  the  main-beam 
direction  8o  is,  from  (A), 

0  <  6o  <  sin”1 (X/2Md)  .  (10) 
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For  «,  and  (10)  glves  0  <  eo  «  2.56-,  a  very  snall  range  ^ 

represents  a  useful  parameter  for  specifying  the  radiation  pattern 
of  arrays  employing  digital  phase  shifters. 


III.  THE  MINIMUM-PHASE-ERROR  CRITERION 


When  the  phase-index  function  f(a,k)  *  0  in  (3),  the  radiation 
pattern  is  identical  to  that  produced  by  continuous  phase  shifters  and 
is  given  by 

K 

Mu>  "  1  \  exp(jku)  .  (11) 

k**-K 


The  difference  between  E(u)  in  (2)  and  E^(u)  in  (11)  represents  the  error  in 
radiation  pattern  due  to  quantized  phase  errors.  If  2Ttf(a,k)/M  <<  1,  the 
factor  exp  [ j  2irf  (a ,k) /M]  in  (2)  can  be  replaced  approximately  by 

I  j27rf (a,k)/M  and  the  mean— square  error  of  the  radiation  pattern  can 
be  written  as 


Eo(u) 


du 


I  (f  (a,k)  }2 

k“-K  K 


(12) 


Equation  (12)  indicates  that  it  is  necessary  to  minimize  the  absolute 
values  of  the  phase-index  function  f(a,k)  in  order  to  minimize  the  mean- 
square  error  of  the  pattern. 

As  illustrations  we  choose  to  compute  for  a  25-element  array  (K=12) 
the  fm(a,k)  distributions  which  have  minimum  absolute  phase-index  values 
xor  two  main-beam  directions  corresponding  to  a  »  1/3  and  a  =  1/2.  The 
results  are  shown  in  Fig.  2.  These  distributions  have  been  used  in  con¬ 
junction  with  the  amplitude  excitations  necessary  to  produce  a  Chebyshev 
pattern  with  a  -40  dB  sidelobe  level  (if  there  are  no  phase  errors). 
Radiation  patterns  have  been  computed  for  M»24»16  and  are  drawn  as  the 
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solid  curves  in  Flgs.  3  and  4.  fThe  significance  of  the  dashed  curves  will 

discussed  in  a  later  section.)  It  is  noted  from  Flg.  3  that,  while  most 

of  the  sidelobes  confers,  to  the  designed  level,  high  sidelobes  appear  in 

two  regions,  the  highest  one  being  only  20  dB  down  from  the  main  beam  at 

-1/3.  The  highest  sidelobe  of  the  pattern  for  a-1,2  in  Flg.  4  i8  _23  dB> 

and  the  solid  curve  hardly  resembles  a  ruac  1. 

y  resembles  a  Chebyshev  pattern.  It  is  evident  that 

the  minimum-absolute-phase-error  a 

P  error  criterion  does  not  yield  very  good  results. 

The  directions  u  where  high  sidelobes  appear  do  not  depend  on  the 
value  of  M,  but  the  levels  of  these  high  sidelobes  vary  approximately 
aversely  wi  h  H.  An  examination  of  the  minimum  phase-index  distribution, 
V°.>0,  in  Fig.  2  reveals  the  presence  of  regularity  in  element  phase 
errors.  This  regularity  can  be  correlated  with  both  the  location  and  the 
level  of  the  high  sidelobes.  A  quantum  (one-step)  change  in  the  digital 
Phase  shifter  of  one  or  several  of  the  array  elements  breaks  this  regularity, 
but  it  does  not  guarantee  a  sidelobe  reduction.  A  modification  of  f  (a,k) 
can  be  found  by  trial-and-error  which  will  result  in  some  sidelobe  reduc¬ 
tion,  but  a  t  ial-and-error  approach  is  impractical  for  an  array  with  a 
large  number  of  elements.  In  the  following  section,  a  systematic  method 
will  be  formulated  for  peak-sidelobe  reduction  by  taking  advantage  of  the 
nonreciprocal  nature  of  latching  ferrite  phase  shifters. 


IV.  NONRECIPROCAL  PHASED  ARRAY 


The  latching  forrite  phase  shifters  useful  in  practical  applications 
are  inherently  nonreciprocal  [3],  and  it  is  possible  to  control  the  phase 
distribution  in  a  radar  array  differently^  for  transmitting  and  for  receiving. 
Hence  sidelobe  reduction  can  be  effected  by  suitably  adjusting  the  phase 
shifts  in  the  elements  of  an  array  for  transmitting  to  counteract  the  high 
sidelobes  in  receiving  through  the  introduction  of  an  appropriate  weighting 
function.  Let  E^(u)  and  E^_(u)  denote  the  transmitting  and  the  receiving 
pattern  functions  respectively,  as  follows: 

Efc(u)  -  l  Ifc  exp(jku)  exp [j  ^  f(a,k)]  (13) 

lc 

Er(u)  -  l  Ifc  erp(jku)  exp[j  ^J-g(a,k)]  ,  (14) 

lc 

where  g(a,k)  is  a  phase- index  function  defined  in  a  manner  similar  to 
f(a,k)  in  (5): 

g(a,k)  «  ak  -  n' (k)  ,  (15) 

and  -2irn *  (k) /M  is  the  digital  phase  shift  of  the  kth  element  for  the 
receiving  state.  When  g(a,k)  »  f^(a,k)  the  mean-square  error  of  the 
receiving  pattern,  Er(u)  -  Eq(u),  is  a  minimum.  However,  high  sidelobes 
appear  in  certain  directions  as  pointed  out  in  the  preceding  section. 

This  drawback  can  be  overcome  by  a  reduction  of  the  sidelobes  in  the 
same  directions  for  the  transmitting  state. 

In  order  to  counteract  the  high  sidelobes  in  the  receiving  pattern, 
we  consider  the  following  mean-square  error  of  the  transmitting  pattern: 


8 


(16) 


2tt 

^  *  2tT  lEt^u)  “  Eq(u)|2  w (u)du  , 

0 

where  w(u)  is  a  nonnegative  weighting  function  which  is  to  have  a  large 
value  in  regions  where  the  receiving  pattern  has  high  sidelobes.  For 
symmetrical  amplitude  and  antisymmetrical  phase  distributions,  (16)  can 
be  expressed  as 


with 


and 


2  K  K 

(TT)  J  £  a,  a  v, 

K  k-i  pii  k  p  k» 


\ '  v(“*k) 


'kp  “  -it  J  sin  ku  sin  Pu  w(u)dti  . 
0 


(17) 

(18) 

(19) 


It  has  been  assumed  in  (19)  that  f(a,0)  -  0,  all  phase  thifts  having  been 
measured  from  the  reference  phase  at  the  center  element.  We  now  consider  the 
minimization  of  the  mean-square  error  £  by  a  one-step  phase  shift  in  the  kth 
digital  phase  shifter  away  from  fffl(a,k),  the  state  where  the  absolute  values 
of  the  phase-index  function  is  a  minimum;  that  is,  we  make 


®k  “  W°’k)  +  ek]  »  (20) 

where 

(  +1  for  f  (a,k)  <  0 

Ek  “  |  (21) 

l  "I  for  fm(a,k)  >  0  . 

When  the  kth  element  is  not  under  phase  compensation,  -  1;  the  phase 
shift  in  the  transmitting  state  is  the  same  as  that  in  the  i 'ceiving  state 
Substitution  of  (20)  in  (17)  yields 
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(22) 


where 


K  K 

[  l  l  VV«IVI^fm^a‘k^fni^a*P^  +  R}  * 

,,  L  ,  kp  k  p  m  m 

k=-l  p=l 


K  K 

7  (V  +  y  v,  e,  )e 
L,  p  . L,  kp  k  p 

p=l  r  k=l 


K 


V  =2  y  v,  I.  f  (a,k) 
P  kti  kP  k  m 


(23) 


(24) 


For  a  given  main-beam  direction  a,  the  first  term  in  the  right-hand 
side  of  (22)  is  determined.  The  sidelobe-reduction  procedure  then  con¬ 
sists  of  searching  for  the  one-step  phase  adjustments  in  particular 
elements  in  order  to  minimize  R,  as  expressed  in  (23) . 


¥.1 _ PHASE  COMPENSATION  AND  WEIGHTING  FUNCTION 


If  a  one-step  phase  adjustment  is  to  be  made  in  the  kth  element 
onJv  in  accordance  with  (20),  we  have,  from  (23), 


R  ‘  \  "  ekVk  +  vkk  •  (25) 

It  is  a  very  simple  matter  to  find  the  value  of  k  which  gives  a  minimum 

R  for  a  specified  weighting  function  by  numerical  computation  of  A  for 

k“l,2, . . . ,K.  For  one-step  phase  adjustments  in  2N  elements,  R  in  (23) 
can  be  written  as 


R  "  l  l  +  2vkpEkep}  ’  <26> 

where  the  summations  for  k  and  p  range  over  the  N  different  integers 
denoting  the  elements  whose  phase  shifters  are  to  be  adjusted.  The 


combination  of  the  N  integers  which  make  R  a  minimum  must  be  found 
numerically,  but  the  computation  is  simple  and  the  procedure  straight¬ 
forward;  there  is  no  need  to  compute  any  radiation  patterns  until  the 
particular  elements  for  phase  compensation  have  been  selected. 

The  choice  of  the  weighting  function  w(u)  for  the  computation 
of  vkp  in  is  flexible.  The  simplest  choice  is  perhaps  to  make 

w(u>  zero  everywhere  except  in  the  regions  where  high  sidelobes  appear 
in  the  receiving  pattern.  We  assume 


w(u) 


1,  u^-Au<u<u^+Au  and  t^-AufuM^+Au 
0,  elsewhere  . 


(27) 
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In  (27)  ux  and  u2  denote  the  locations  and  Au  the  extent  of  the  sidelobea 
to  be  reduced.  Examination  of  the  solid  patterns  in  Figs.  3  and  4  reveals 

that  it  is  appropriate  to  choose  u2  =  2,^.  With  this  choice,  we  obtain, 
in  substituting  (27)  in  (19), 


and 


cos (k-p)u^ 


sin (k+p) Au 
(k+p) Au 


cos (k+p)u^] 


(28) 


sin  2kAu 
2kAu  cos 


2k  u,  ]  . 


(29) 


The  above  expressions 


are  convenient  to  use  in  computing  R. 


VI.  ILLUSTRATIVE  EXAMPLES 


The  weighting  function  in  (27)  has  been  used  to  reduce  the  high 
sldelobes  which  appear  in  the  radiation  pattern,  Er(u),  of  digitally 
controlled  arrays  designed  under  the  minimum-absolute-phase-error  cri¬ 
terion  in  Section  III.  From  Figs.  3  and  4  it  is  ~een  that  U;1  -  2tt/3 
and  9 it/10  respectively  for  o  -  1/3  and  1/2.  An  appropriate  choice  for 
is  tt/10  for  both  cases.  The  essential  data  are  presented  in  Table  I. 


TABLE  I 

Phase  Compensation  Data  for  25-Element  Array 


a 

1/3 

1/2 

w(u) 

ul  "  2ir/3,  u2  -  Att/3 

Au  -  tt/10 

ux  -  9 tt/10,  u2  -  llu/10 

Au  «  tt/10  J 

Case 

(A) 

asa 

(E) 

(F) 

(G) 

Pai^s  of  elements 
compensated  (N) 

1 

2 

3 

D 

2 

2 

2 

Eleme 

nt  Nos.  (k) 

1 

wm 

1,8,10 

5 

m 

3,9 

1,5 

R 

Value 

-0.3301 

min. 

-0.3783 

min. 

L 

-0.3791 

min. 

-0.5503 

min. 

-0.6549 

min. 

-0.6193 

2nd. 

-0.6125 

3rd.  | 

Rank 

In  case  (A)  for  a=l/3,  only  one  pair  of  elements  (N=l)  are  phase- 
compensated,  and  calculations  with  (25)  show  that  R  is  minimum  when  k-1. 
irom  Fig.  2(a),  it  follows  that  the  phase  in  the  first  element  to  the 
tight  of  the  array  center  should  be  shifted  down  by  one  step  and  the  phase 
in  the  first  element  to  the  left  of  the  array  center  (k-1)  should  be 
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shifted  up  by  one  step  on  account  of  the  assumed  phase  antisymmetry.  The 
equivalent  pattern  defined  by 

E(u)  =  |Er(u)  Et(u)|1/2  (30) 

for  case  (A)  is  plotted  as  the  solid  curve  in  Fig.  5.  In  comparison  with 
the  solid  curve  in  Fig.  3,  it  is  seen  that,  although  reduced  in  level  by 
slightly  more  than  A  dB,  the  high  sidelobes  remain  visible  at  about  the 
same  locations. 

If  two  pairs  of  elements  are  compensated  (N=2)  ,  the  phases  in 

elements  k*=A  and  k=5  should  be  adjusted  one  step  down  and  up  respectively 

from  the  f  (1/3, k)  values  shown  in  Fig.  2(a)  in  order  to  obtain  a  minimum 
m 

R.  This  is  listed  as  case  (B)  in  Table  I.  The  resulting  equivalent 

pattern  is  plotted  as  the  dashed  curve  in  Fig.  5,  which  exhibits  a  further 

« 

reduction  in  the  levels  of  the  high  sidelobes.  In  case  (C),  one-step  phase 
adjustments  are  applied  to  three  pairs  of  elements  (N=3) .  The  excitation 
phases  in  the  1st,  8th,  and  10th  elements  are  shifted  down,  up,  and  down 
by  one  step  respectively  for  a  minimum  R.  The  equivalent  pattern  is  shown 
as  the  dashed  curve  in  Fig.  3.  Comparing  with  the  solid  curve  represent¬ 
ing  the  receiving  pattern  in  the  same  figure,  we  see  that  the  highest 
sidelobe  has  been  reduced  from  -20  dB  to  -29  dB,  an  improvement  of  9  dB, and 
that  the  low  sidelobes  in  other  regions  have  been  raised  slightly. 

A  similar  procedure  is  applied  to  reduce  the  high  sidelobes  in  the 
receiving  pattern  for  a=l/2  in  Fig.  A.  Four  cases  are  studied  using  a 
weighting  function  w(u)  as  specified  in  (27)  and  Table  I.  In  case  (D) , 

N«1  and  k®5,  the  phase  in  the  5th  element  is  shifted  down  by  one  step  for 
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nlnin,™  R.  Thc  equivalent  pattern,  plotted  as  the  solid  curve  in  Flg.  6, 
still  shows  the  existence  of  high  sidelobes  on  both  sides  of  u-r.  Case 
for  phase  compensation  in  tvo.pairs  of  elements  (N=2,  indicates  that  one. 
step  phase  adjustments  are  required  in  the  3rd  and  7th  elements  in  order  to 
fake  R  minimum.  The  corresponding  equivalent  pattern  shown  dashed  in  Fig.  6, 
however,  reveals  that  the  original  high  sidelohes  have  been  over-reduced 
At  the  same  time,  the  sidelohes  in  some  other  regions  are  made  unduly  high 
m  an  effort  to  correct  this  situation,  the  cognations  of  4-eleme„t  (N-2) 
Phase  compensation  to  give  the  second  and  third  smallest  R  „ere  examined  '  ' 
These  are  listed  as  cases  (F)  and  (G)  in  Table  1.  Pattern  calculatIons 

indicate  that  case  <F) .  like  ease  (E) ,  still  over-reduces  the  high  sidelobes 

near  u-s  at  the  expense  of  raising  the  sidelobes  in  other  regions.  For  this 

reason  the  pattern  for  case  (F)  is  not  shown.  Case  (C)  with  k-1  and  5  proves 

to  be  a  good  choice,  aiding  the  dashed  pattern  in  Fig.  4.  Comparison  o, 

the  two  patterns  in  Fig.  4  discloses  a  reduction  of  5  dB  in  the  p,ak-sidelobe 
level. 

Obviously  the  same  technique  can  be  extended  to  achieve  further  re¬ 
duction  in  sidelobes  b;r  phase  compensation  in  more  than  two  pairs  of  ele- 

Ultimately  it  is  a  compromise  between  the  amount  of  further  improve- 
merit  and  the  attendant  complexity. 
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VII.  CONCLUSION 


It  has  been  shown  that  the  radiation  pattern  of  a  digitally  con¬ 
trolled  phased  array  can  be  cade  periodic  with  respect  to  a  now  parameter. 
In  so  doing,  the  study  of  such  an  array  is  much  simplified  and  it  is  only 
necessary  to  consider  the  pattern  within  a  very  narrow  range  of  the  scan 
angle.  The  simple  expediency  of  minimizing  the  mean-square  error  of  the 
pattern  function  does  not  yield  uniformly  low  sidelobes  throughout  the 
visible  legion.  By  taking  advantage  of  'he  nonreciprocal  nature  of 
latching  ferrite  phase  shifters,  a  systematic,  easy-to-apply  method  has 
been  developed  for  the  reduction  of  ponk  sidelobes.  The  same  technique 
can  be  used  to  suppress  interference  or  improve  signal-to-noise  ratio  in 
radar  applications.  All  that  is  necessary  is  to  choose  a  weighting  func¬ 
tion  w(u)  which  has  a  large  value  in  the  regions  (directions)  of  inter- 
ference  or  high  noise. 
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MAIN-BEAM 

DIRECTION 


Fig.  1  -  A  linear  array  with  uniform  spacing. 


fm(a,k) 
I 

2/3 
1/3 
O' 
-1/3 
-2/3 
-I 


T~'  — ' 

4  5  6 

(a) 


8  9  10  II  12 


* 


Fig.  2  -  Minimum  phase-index  distributions  for  a  25-element  array 
(a)  a  =  1/3,  (b)  a  =  1/2. 
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